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On the Deflection of the Vertical by Tidal Loading of the 

Earttis Surface. 

By Horace Lamb, F.R.S. 
(Eeceived January 25, 1917.) 

The earliest calculations on the above subject in which the flexure of the 
earth's surface is taken into account, as well as the direct attraction of the 
water, are contained in a well known pager by G-. H. Darwin.* It was 
there shown, in a particular case, that the tilting effect would conspire with 
the attraction, and would on certain assumptions be directly proportional to 
it. The latter remark was stated as due originally (in a more general form) 
to Sir W. Thomson. On inserting numerical values it was found that the 
apparent deflection due to tilting would considerably exceed that due to 
..attraction.")- The whole subject has of late excited renewed attention, owing 
to its bearing on observations of lunar deflection of gravity ,| and in a 
recent paper by Terazawa§ the matter is specially considered from the latter 
jpoint of view. 

The present paper, after discussing a few typical problems, goes on (in § 4) 
to examine the effect of one or two considerations which have hitherto been 
disregarded in such calculations. It is true that the corrections involved 
are under some conditions negligible, but they are of theoretic interest, and 
it is found that at great distances from a load, an4 therefore in all cases of 
widely distributed load, they may attain considerable relative importance. 
In the first place, owing to the deformation of the earth's surface and the 
altered distribution of density an additional horizontal component of force 
on the plumb-line is introduced ; this tends to counteract the attraction of 
the water. A more important point is that the influence of gravity on the 
deformation has been ignored. In attempting to estimate the effect of 
gravity it has been found convenient, in order to avoid difficulties not 
altogether of a mathematical kind, to limit the investigation to the case of 
incompressibility. I have also neglected, in the first instance, the dis- 
turbance in the field of gravity itself, due to the load and the deformation, 
so far as this 'affects the strains. When subsequently the alteration of the 

* ' Brit. Ass. Keports,' 1882 ; < Scientific Papers,' vol. 1, p. 444. 

t The relation was pointed out independently by Chree, 'Phil. Mag.' (5), vol. 43, 
p. 177 (1896). 

% Love, 'Some Problems of Geodynamies,' p. 88 (1911). 

§ " On Periodic Disturbance of Level Arising from the Load of Neighbouring Oceanic 
Tides," ' Phil. Trans.,' A, vol. 217, p. 35 (1916). 
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field is taken into account a curious point arises. Eor mathematical 
simplicity the " earth " has been treated, as is usual in such investigations, 
as flat and infinitely extended. It appears that if this were the case the 
surface would be unstable, whatever the degree of rigidity, for disturbances 
exceeding a certain critical wave-length. This wave-length is, however,, 
enormous, and reason is given for the view that inferences can still 
legitimately be drawn from our results as to the character of the effects 
actually produced. 

1. For the present purpose, the earth is regarded as a semi-infinite 
isotropic elastic solid, extending indefinitely on the lower, positive, side of 
the plane 3 = 0. If a load of surface density 

a- = Jo (mr), (1) 

where r, = ^/O^ 2 + 2/ 2 )> denotes distance from the origin, be placed upon it, 
the usual theory, in which elasticity only is taken into account, gives a 
surface depression 

Wo = &=lHl . J o (""P. (2) 

fi m 

and a dilatation 

A = - iSkl^Ill e ~™ J (mr), (3) 

where fi denotes the rigidity, and tj Poisson's ratio.* 
The gravitation potential due to the load (1) is 

V = 2 ^Jo(^) ; (4 ) 

m 

where 7 is the gravitation constant, and the upper or lower sign in the index 
is to be taken according as z is positive or negative. For the two values of 
V satisfy Laplace's equation v 2 V = ; they are continuous at the surface ; 
and they give the proper discontinuity in the value of dV/dz. The surface 

potential is therefore 

t 7 o Jo (mr) /KN 

Vo = 27T7-^ — '-. (5) 

m 

The apparent deflection of the plumb-line due to tilting is —dwo/dr, and 
the deflection due to the attraction of the load is —dV /gdx. The ratio of 
the latter to the former is 

2tt7//, _ 3 p/p 

(1-V)g 2 2(1-10' gR' 

where E denotes the radius, and p the mean density, of the earth. Since 
this ratio is constant and independent of m, we infer that it is the same for 

* 'Proc. Lond. Math. Soc.,' vol. 34, p. 278 (1902). The formulae are slightly simplified 
by the use of rj y = iV(^ + / x )> in place of X. 



(6) 
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all symmetrical distributions of load, and therefore for a concentrated load, 
.and therefore for all distributions whatever. This verifies Thomson's 
remark already referred to. If we take for numerical illustration* 

97=-£, ^ = 4xlO n , = 981, 7=6*66x10-8, 
in C.G.S. units, we find 

jp^ . = 0-232. (7) 

In the case of incompressibility (?? = |) we have, with the above value of /i, 

^ = 0-348. (8) 

g 2 

The method of passing from (1) to the case of any symmetrical distribution 
of load is explained in the paper referred to. To deduce the effect of a con- 
centrated load M at the origin we multiply the various expressions by 
M.mdm J 2tt y and integrate from m = to m = co . Since, however, the 
deflection due to the attraction of M is jM/gr 2 we have at once 

_dW0 = (I"" q)ff ^ M ,gv 

dr 2tt{j, * r 2> 

in virtue of the proportionality referred to. 

2. In a case specially discussed by Terazawa the load is supposed to be 
distributed uniformly over a circular area. This can be treated very simply 
with the help of a formula given by Ferrers, f He proves (in effect) that the 
potential of a circular disc of radius a and surface density 



/yta\ 7b 



"= !-^ > ( 10 ) 



a 



&t external points in its plane is 

n = m^jV a? 



f* qq / __^ \ *n *4- Jt 



r 2 \ n+ % d\ 



r (^.+l) J V ^ 2 +V (^ 2 +x)v / x , 

,the lower limit of integration being r 2 —a 2 . 

In the present case we have n = O.J Putting a 2 + \ = r 2 /sin 2 <9, we find 

°» = 4a2 tv¥^^e) = T ^ " (1 - /fi2 > *>• ^ 

* It is very uncertain what are the best representative values to take for /* and p 
•(the density), on account of the actual variation with depth. If we put p = 2*5 the 
above value of p makes the velocity of propagation of transverse vibrations agree 
with the value found by seismologists for the upper strata, viz., 4 km. per second. 

f " On the Potentials of Ellipsoids ... of Variable Densities," i Quart. Journ. Math.,' 
vol. 14, p. 1. 

J In the case n = J, also treated in a different manner by Terazawa, the integral 
reduces ; thus 

Q o = «« { (i-| 5) sin " 1 f+ r v (5- 1 ) }• 
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in the notation of elliptic integrals, the modulus being 

k = a/r. (13> 

For the attraction we derive the simple expression 

- ^ = ^Fx-Ex). (14> 

Tables of Fi and Ei, abridged from Legendre's, are included in several modern 
collections. 

The object of Terazawa's paper was to form some estimate of the possible 
effect of the periodic tidal loading and unloading of the Atlantic area on 
Michelson's observations of lunar deflection of gravity at Chicago. The radius 
of the circle was taken as 2000 km., and the distance of the point at which 
the effect is to be estimated as 3000 km. If we put k = | we find from 
(14) on reference to the Tables 

dr 

To find the deflection of the pendulum by attraction we must multiply this 
by ypihfg, where pi is the density of the water, and h the tidal elevation,, 
Assuming p x = 1, h = 100, we find 

- i^l = m x 10 ~8 ; ( 16 ) } 

gdr 

which is Terazawa's result, obtained by a different process. The tilting effect,, 
obtained by multiplying by 4*31, is 5*04 x 10" 8 . 

Since tidal elevation in one part of the ocean is compensated by depression 
elsewhere, it may be worth while to notice the case where the tidal elevation' 
in a circular basin is supposed to vary as the distance from a nodal line 
(x = 0). It follows from another formula of Ferrers that the potential, at. 
points in its plane, of a disc whose surface density is x/a, is 

XV = 2 A f(l-^) \ a2+ ^ Vx > (17> 

the lower limit being r 2 — a 2 . If we put x = r cos <j>, this is equivalent to 

4 a 3 . ft" sin 2 # cos 2 # d6 



f*" 

♦f. 



Qo' - r ^T^^y^^^^^^y 



l^i{(2-F)E 1 -2(l-F)F 1 } > (18) 



3 Je 2 
where k = a/r as before. From this we find 



3^ = 4c^ {4(Fi _ Ei) _ F(Ei + Fi)} _ (19) 

cr 6 k 
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If the height of the water be hx/a, the radial deviation of the plumb-line 
due to attraction is found by multiplication by ypih/g* If pi = 1, A = 100, 
Je = -§-, I find 



av, 



__ Q.11 n ~ a A. ^ 1 A-9 



gor 



3-11 cos ^xlO" 9 . (20) 



(22) 



The tilting effect is therefore, on the previous assumptions, 

_ dwo _ 1 .g 4 cos . x 1Q - 8> . 2 1) 

or 

3. Two-dimensional problems of the above kinds can be treated ah]initio by 
a simpler analysis, but we may conveniently use the result contained in (9), 
The loads are now supposed applied on infinitely long bands of the surface, 
parallel (say) to the axis of y. Unless the total load (positive'^and negative) 
is zero, the integrals which give the surface depression and the potential are 
divergent ; but we may still draw inferences as to the deflections produced by 
bands of great length at distances small compared with the length. 

Thus, for a linear load of amount M' per unit length, we have 

_dY _ 2 7 M / 

(JCVJU JU 

by the theory of the attraction of an infinite straight line, and therefore 

dw _ (l — y)g M^ ^\ 

dx 7T/JL X ' 

If we imagine a band of breadth 2a to carry a uniform load pji per unit 
area, we have f or x >a 

_ d Jl = 2 VP lh log ^±± , (24) 

gdx g x—a 

whence the value of — dw /dx follows as before. 

A more interesting case is where the load varies as the distance, positive or 
negative, from the medial line x = 0. For the sake of comparison with a 
problem discussed by Darwin,* the depression within the band, as well as at 
external points, is considered. The load per unit area being pihx/a, for 
x 2 < a 2 , we find from (23) 

^ w o _ (1— ^)#piA f a x'dx' 
dx irfia 

If x > a this gives 






(25) 



dwo _ (1— rj)yp x li lx lo „£+a_ 2 \ 



(26) 



dx TTfj, \a a x—a 

For 0<«<awe have, taking the " principal value " of the integral in (25), 

dw Q _ (l—7))gpih /x lo „ a+x _o\ / 2 1) 

dx tt/jl \a a a— x / ^ * 

* Loc. cit. ante. 
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Hence, for x > a, 

(l—r))gpihafx x^—a 2 , x-\-a\ /oox 

W ° = H \a—2^r l0 Z^ (28) 

and, for < x < a, 

= (l-v)gi P haix + M 1 a+xy 
TTfi \a 2a J a— xj 

The arbitrary constants of integration have been adjusted so as to make w 
vanish for ^ ~ oo , and to establish continuity at the edge x = a. Since iv 
must be an odd function, the above expressions will hold also for % negative. 
The graph of w Q is shown in the figure. 



The tilting effect at external points is given by (26) ; from this the direct 
attraction of the load can be deduced. 

4. Besides the direct attraction of the load on the pendulum there will be 

a horizontal force due to the changed configuration of the solid earth. If 

we assume the formulae (2) and (3) for the deformation, the corresponding 

potential is found to be 

_ HE . J o( my ') ? (30) 

fi m 2 

which gives in the case of a concentrated load a deflection 

-^ M. (31) 

2 /uL r 

The effect on the plumb-line is therefore equivalent to a repulsion varying as 
the inverse distance, whereas the attraction of the load varies as the inverse 
square. It is not worth while to give the proof of the formulae, for the effect 
could only become appreciable at very great distances from the load, and in 
these regions the deformation, which has been assumed to be given by (2), is 
seriously affected by the action of gravity, which is next to be considered. 
The effect here referred to is included implicitly in the amended results. 

5. When we proceed to examine the effect of gravity in modifying the 
deformations, we are met by difficulties relating to the state of initial stress 
of the earth. To avoid these the material is here assumed to be incompres- 
sible, and the initial stress to consist of a uniform hydrostatic pressure* 

Pi = gp%. (32) 

* Of. Rayleigh, "On the Dilatational Stability of the Earth," 'Roy. Soe. Proe., 5 A, 
vol. "77, p. 486 (1906) ; < Scientific Papers, 5 vol. 5, p. 300. 
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It is also assumed for the present that in calculating the strains we may 
neglect the alteration in the distribution of gravity due to them and to the 
load. 

It is convenient to begin with the two-dimensional case. On the above 
hypothesis the elastic equations are, in the. usual notation, 



fdho 8% 
\dx 2 dz 2 , 

fc 2 w , 8V 
\dx 2 dz 2 



dp 

dx 

dp 

dz 



^N 



>> 



J 



where p denotes the excess of pressure over its initial value p\. 
Assuming for a typical solution, 

p = Ce~~ mz cos mx, 
we have 



% = (z + A)e~ mz sin mx 



~\ 



to = — + B)<r 

The condition of incompressibility, viz, 

du s dw _ 
dx dz 



y 



•mz 



cos mx 



j 



gives 



A-B = 



0, 

_1 



whilst the condition of vanishing tangential stress, viz. 



du dw __ n 
dz dx 



at the surface z = 0, gives 



A + B = 



m 



Hence A = 0, B = Ijm, so that 

C 



~\ 



u =. 



mz e mz sm mx 



■ID 



fjum 

C 



>• 



(l + mz)e mz cosmx 
2/nm ^ 



It may be noticed that divfdz vanishes at the surface. 

If a be the surface-density of the load, the surface condition is 

dw 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



■Pi-P + 2/j, 



iz 



■ga, 



(41) 



VOL. XCIII. — A. 
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to be satisfied for z = iv . Since dpijdz = #/>, this becomes 

#P^o +p — 2 ^ — = ga } (42) 

to be satisfied for z = 0. Hence if 

a- = H cos ma? (43) 

Tims Wo = §5 . 22i^, (45) 

where # = f/p/2/x. (46) 

To deduce the effect of a concentrated load of line-density M' we put 
H = M / dm/ir, and integrate from to oo . Thus 



we have 



M'# f cos mo? 7 M# f 00 me ro;c 7 //i)7 v 

w = _*- . — dm = 7—^- -tt: — s rfm, (47) 



* 



27r/xJ m + /c 2 777* Jo m 2 ~f /c 2 



or w == 77— { ( i ^ ~" Si «) sin /&c — Ci kx cos ##} , (48) 

in the ordinary notation of the sine- and cosine-integrals, viz. 



Six 



rx smx -, „. f 00 cos# 



rf#, Qix ~ — I cfa?. (49) 



^ J a: "^ 



We may notice that the value of w is now finite except for x = 0. From 
(48) we deduce 

-^:=^A-K(i7r-Si^ (50) 

aa? 27rjC6 {jo J 

which now takes the place of (23). 

As regards the numerical value of k we have, taking //,/p = 16 x 10 1Q , 

k~ 1 = 3260 km. (51) 

In any practical application we are therefore concerned only with moderate 

values of kx. 

The following table gives in Column I values of the function 

— - dm = (W~Si/caAeos/c^-f Ci/c$eos/a£ ? (52) 

based on Glaisher's tablesf of the standard integrals (49). In Column II 

* As to the equality of the two integrals, see the author's ' Hydrodynamics ' (1916), 
Art. 243, where other references are given. 

t * Phil. Trans.,' p. 367 (1870). An abstract is given in Dale's 4 Five-Figure Tables 
of Mathematical Functions,' London, 1903. 
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this function is multiplied by koc, and the result subtracted from unity. This 
shows the ratio in which the tilting effect is diminished from its value as 
given by the formula (23), where gravity was neglected. At small distances 
the change is slight, but at a distance of 1000 km., for example, the tilting 
is reduced by the action of gravity to the extent of 30 per cent. 



fC*U* 


I. 


n. 


ft%&* 


I. 


II. 


o-o 


00 


1-0 


0*5 


-8605 


0-570 


0*05 


1 -3980 


0*930 


0-6 


-7985 


-521 


o-i 


1 -2910 


0-871 


0*7 


-7451 


0-478 


0-2 


1 -1369 


0*773 


0*8 


*6987 


0-441 


0-3 


1 -0236 


0*693 


0-9 


'6578 


0*408 


0-4 


'9341 


0*626 


1*0 


*6214 


0-379 



It is to be remarked that on the present suppositions the residual 
horizontal attraction on the pendulum, when allowance is made for the 
altered attraction of the solid, bears the same constant ratio to the tilting 
effect as in Art. 1 above. Eecurring to the distribution (43), and assuming for 
the potential due to the load and to the change of configuration of the solid 



V = 'De Tmz cos mx, 
we have, on account of the discontinuity in SV/dz, 



(53) 



2mJ) cos mx = 47ry (<r~ pw^) y 



whence 



Thus 






27r*yH 

m + tc 



Y _ Airy/u 



gw 



9' 



(54) 
(55) 



(56) 



Since this does not involve m, the ratio will hold for all two-dimensional 
distributions of load. With the former numerical data its value has been 
found to be 0'348. 

6. A similar investigation can be made for the case of a load concentrated 
at a point. The elastic equations are 



s dp 



fM^H 



dp 
df 



fl^ 2 W 



dz 1 



involving also 



VT = 0. 



(57) 
(58) 



If q denote radial displacement in a horizontal plane, so that 

n = qxjr f v = qy/r, 



(59) 



Jj A A 
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the equations become, in the case of symmetry about the origin, 



\dr 2 r dr r 2 dz 2 / dr 



M 



fd 2 w ldw d 2 w\ ___ dp 



s 






lor 2 ' r dr ' 3« 



+ 



with 

Assuminj 
we find 



d 2 p 1 dp d 2 p 

dr 2 t dr dz 2 



"&2 1 dz 
= 0. 



^> = Ce WMr J (mr), 

C 
([ = -£- (z + A) e~ ^ Ji (mr) 



,/x 



w 



2yU, 



y- 



(z + B) g mz Jo (mr) 



.y 



The condition of incompressibility, viz. 

dq q dw _ 

r r dz 

A—B = — 



0, 



leads to 



1 



m 



whilst the condition of zero tangential stress at the surface, viz. 

1 



dq _Jdw 
dz or 



gives 



A + B — 



m 



Hence A = 0, B = 1/m, so that 

C 



■ID 



2/JL7K 



(1 -f mz) e~ mz Jo(mr). 



The surface-condition has again the form (42). If the load be 

a = HJo(mr) 



we find 



w = 



gH Jo (mr) 



(60) 



(61) 
(62) 



(63) 



(64) 
(65) 



(66) 
(67) 



(68) 



(69) 

(70) 



2jju m + tc 
where k has the value (46). 

To pass to the case of a concentrated load M at the origin, we put 
H = Mm<:Zm/27r, and integrate from to oo . Thus 

Mg n 



47774 J m + fc 47TyC6 [_r 



»,*^<to} 



! m-\- tc 

Mg j 1 



■4 



4717-1 Lr 



«: 



00 ] 



* A proof of the equality of the two integrals last written is given in the Appendix 
to this paper. 
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An asymptotic expansion of the latter integral is easily obtained, but has 
hardly any application from the present point of view. For moderate values 
of kt we have recourse to the formula* 



poo 

l e ' 



■xsmhu 



du 



•iir 



sin (x cos 6) d6 + \ it K (x), 



(72) 



JO 

where K (x) is the BesseFs function " of the second kind "denned by 



l 00 



K (x) = — cos (x cosh u) du. 

IT J 



(73) 



Hence 



Jo(mr) 7 (/cr) 2 

- - n.rn. = *-"}* — — 



+ 



(/C^)' 



>-...+ A 7rK (*r). (74)f 

The function j7rK has been tabulated by B. A. SmithJ under the designa- 
tion /eJo — Y , where k stands for log 2 — C or 011593, C being Euler's 
constant. With this assistance the following Table has been constructed. 
Column I gives the values of the function (74); in Column II these are 
multiplied by kt and subtracted from unity. The result gives the ratio in 
which the surface depression w is diminished by gravity : — 



Kr. 


I. 


II. 

I 


Kr. 


I. 


II. 


o-o 


GO 


i 
1-0 i 


0-5 


1 -1845 


0-408 


0*05 


3 '1591 


0-842 


0-6 


1 -0609 


0*364 


o-i 


2 -5098 


0-749 


0-7 


-9621 


0-327 


0*2 


1 -8973 


0-620 


0'8 


-8807 


0-295 


0-3 


1 -5650 


0-530 


0-9 


-8125 


0-269 


0-4 


1 '3448 


-462 


1-0 


-7541 


0-246 



The tilting effect is 



dwo M# J 1 



dr 



^ir/jb 



K 



rcc m Ji (mr) 
o m + te 



dm 



(ft 



From (74) we have 

1 f °° m Ji (mr) , 
kJo m + ic 



|/cKi(/cr) 



1 2 .3 + 



(kt) 4 



1 2 .3 2 .5 



(76) 



where Ki is written for — K</. The function J77-K1 has also been tabulated 
by Smith under the designation /cJi — Yi. The following Table gives in 

* Of. ' Hydrodynamics,' Art. 194. 

f This formula occurs in Nielsen's ' Handbuch der Cylinderfunktionen,' Leipzig, 
1904, p. 238, but the derivation is different. The ascending series occurs in another 
connection in Kayleigh's ' Theory of Sound,' Art. 302, where it is denoted by j7rK . 
Nielsen denotes it by ^ttZ , and writes — Y for the function on the right-hand side of 
(73). The notation in the text is .employed by H. Weber. 

J ' Phil. Mag.' (5), vol. 45, p. 122 (1898). 
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Column I the values of the function in (76); in Column II these are 
multiplied by (kt) 2 and subtracted from unity. The results show the ratio 
in which the tilting effect is reduced by gravity :— - 



KT, 


I. 


II. 


KT. 


I. 


— ^-^,.....„^. MW , I ,„ I , M . «.„...— .,-,..,1, .nn-n. 

«■ 

JL JL • 


o-o 


00 


1-0 


0-5 


1 '3934 


0*652 


0-05 


19 '0911 


0-952 


0-6 


1 -0971 


-605 


o-i 


9 -1490 


0-908 


0-7 


-8914 


0-563 


0*2 


\ -2342 


0-831 


0-8 


-7418 


0'525 


0-3 


2 -6318 


-763 ' 


0-9 


-6293 


0-490 


0-4 


1 -8501 


-704 ! 


1*0 


-5426 


0-457 



It is easily seen that the residual horizontal attraction will bear to the 
tilting effect the same constant ratio as in the two-dimensional case of the 
preceding section. 

7. It remains to take into account the fact that the gravitational field is 
altered by the presence of the load and the deformation of the surface, and 
that the strains are thereby, to some extent, modified. 

Taking first the two-dimensional case, the equations (33) are replaced by 



ld 2 w , 



d 2 u 



d 2 w 



\ dx 2 dz 2 , 



dx dx 

^dp_ 9V 
dy dz 



s 



(77) 



J 



where V is the potential due to the load and to the change of configuration 
of the surface. If we assume 



p = Ce~ w * cos mx 
as before, and pY = Cie^^cosmx, 

we have, writing C — Ci for C in (40), 

w Q = -— ~ — - cos mx, 

2JJL971 

The load being a = H cos mx, 

the surface-condition (42) gives 

#^(C-Ci) + C = gH. 

2/jum 



(78) 
(79) 



(80) 
(81) 



(82) 



Since the value of dY/dz just above the surface must exceed the value just 
below by 47ty(o-— pwo), we have also 

2j^ = 4^fH-£ig=£ri. (83) 

p L 2/tm J 
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Hence d = 2?r W m C (84) 

and (m 2 + S£-m — ^^- ) C = m%H. (85) 

Hence, putting 

2 

2 /a //, 

we find w = tt— • m ~ a ffM a zH2 cos mx 9 (87) 

2/x, (m-\-cx){m — p) 

V 27T7H m /oox 

-J! = 1 — . ._. __ cos mx. (%$) 

g g (m+ct)(m — /3) 

It appears that for m = /3 the equilibrium would be neutral ; i.e. a 
deformation of wave-length 27r//3 could subsist independently of any- 
external load. For greater wave-lengths the equilibrium would be unstable.* 
This, gravitational instability of an infinitely extended plane surface may be 
unexpected, but it is easily verified in the case of a plane sheet of water, 
where /j, = 0. It is evident that in the absence of the earth's attraction and 
of rigidity the mass would tend to cumulate. The introduction of ordinary 
gravity, and even of rigidity, merely imposes a lower limit to the wave-length 
at which instability sets in. 

This limit is however very high. From (86) we find 






> ■ (89) 



J 



"With the previous numerical values of the constants I find 

a" 1 = 2560 km., /3~ l = 11,940 km., (90) 

whence 27r//3 = 75,000 km., or nearly twice the circumference of the earth! 
This indicates that the paradox has arisen through our ignoring the curvature 
and restricted dimensions of the earth. 

Waiving this difficulty for the moment we have 

-^ = ff{ 1 -^( w ^ + ^)} sin ^ (91 > 

gax g L a + p\m-fa m— -ft/ J 

* This is seen most by directly investigating the effect of a surface pressure P cos msc, 
instead of a material load. We should find 

P m 2 

w = — - . — cos mx. 

which has the same sign as P, or the opposite, according as m < 0. 



306 



Pror. H. Lamb. On the Deflection of the 



Hence, for a concentrated linear load,* 



dw _ M f g f" 1 ac 



dx 4:7Tfjb [_x a 2 — j3 



2 



00 sin mx 
m 4- # 



dm 



£ s 



•qo 



sin ma; 



a 



■yS 2 Jo m— /3 



dm 



(93) 



d Y _ 7M' f 1 a 2 f °° sin mg , 



/3 2 f °° sin m$ ^ 



a + /3J m — 13 J 



dm ^ . (94) 



The second definite integral in these expressions is indeterminate to the 
extent of additive terms of the form 



A cos f3x -h B sin fix. 



Its " principal value " is 

~°° sin mx 



P 



•oo 



o ( m—(3 



dm = it cos fix- 



sin mx 



dm 



o m + /3 

— (j7r4-Si/3^)cos/3^— Gif3x8mf3x. (95)-J- 

The adoption of the principal value of an indeterminate integral may 
appear somewhat arbitrary, and it may he questioned whether results 
obtained in this way have any physical significance. The consideration of 
the following modified problem may perhaps serve to justify the procedure. 

Suppose that we have an endless series of alternately positive and negative 
loads + M' at the points x — 0, + a, ± 2a,.... These may be repre- 
sented as the limit when h-+l of 



M'/? 7TX , 7 o OTTX , 7F; 5tfX 

a = — (h cos h h cos — 

a \ a a 



-f h b cos 



a 



-j- . . . 



(96) 



For, summing the series, we have 

M'A 



l-# 



O" 



a 



l-2A 2 cos — + A 4 



(97) 



the limit of which is zero, except when x is zero or a multiple of ±a. Also, 
if b < a we have 



* b 7 2M' . r27fcsin(7r6/a) 

o-dx — tan x ^ — '— 

-b 7T 



1-h 2 



(98) 



the limit of which is M', however small b is taken. Similarly we should find 



ra+b 

lim. o-dx = — M 7 , 

JU-6 



(99) 



* The term l/# arises as 



See the Appendix. 



lim. 



00 



e-m* sin mxdm. 
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and so on. Now, referring to (91) we have 

dw = MV £ J" x * 3 1 £L_ 1 I ^sin 57 ^, 

dx 2/jua s L ex 2 — fi 2 sir j a -\- a ot 2 —fi 2 srr/a—fi ) a 

(100) 
where s = 1, 3, 5, .... Here 

A (1+ A 3 ) sin — 

2 A* sin— = : k — , (101) 

* l--2A 2 cos — + A 4 

a 

the limit of which is | cosec (irx/a). In the remaining terms we may put 
h = 1 at once. It is shown in the Appendix that for values of # between 
anci «■ a 

v sin (sTTx/a) oca f 00 cosh (J a— a:) 77 , n09\ 

s S7r/cc + ot ~ 27rJo (?7 2 + a 2 ) cosh J a?7 

v sm(s7rx/a) __ 1 cos /3 (| a •—#) /3a P°° cosh (|a — 2;) 7? ^^ 

s sirja—fi ~~ 2 cos -| /3a 27rJ (?? 2 + /3 2 ) cosh |- a?/ 

It is not difficult to see that if /3a be moderately large, but less than it, whilst 
fix is small, the result tends to the form 



dwo _ M.'g 



1 



OL 



3 



.0? a 2 —fi 2 i 



00 fltg-q* /3 3 f f 00 fie-^dr) 

oV 2 + * 2 V * 2 ~fi 2 \_ 7rC08 A^ — J o ^ + /3 2 



(104) 
which is equivalent to that above obtained. 

If we impose periodicity there is in the present problem no indeterminate- 
ness, and no instability so long as /3a < 7r. If we imagine a cylindrical 
earth, periodicity is essential in the solution, and there is no mathematical 
absurdity in regarding the points x = 0, x = 2a, for instance, as being the 
same place. The curvature would, of course, involve some modification of 
the elastic equations. 

I have made some calculations illustrative of the formulae (93), (94). 
To simplify the numerical work I have assumed fi = lot, which is not quite 
consistent with (90), but is equivalent to diminishing the value of /j, by 
about one-tenth. In Column I below the tilting effect (—dw Q /dx) is com- 
pared with that given by the formula (23) where gravity was altogether 
neglected. Column II gives the ratio of the total deflection 

cl 



CvX 



(wo + V /#) 



to the value given by (23). It is remarkable that the numbers in Column II 
only differ by about 1 per cent, from those which we should obtain by 
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ax. 


1-0 


II. 


ax. 


I. 


11. 


o-o 


1-3125 


0-5 


-5503 


0-760 


0-05 


-9265 


1 '222 


0-6 


-4912 


0-691 


O'l 


'8641 


1-145 


0-7 


-4457 


0-636 


0-2 


. -7602 


1-017 


0*8 


*4050 


0-588 


0-3 


-6757 


0-914 


0-9 


-3655 


0-542 


0-4 


-6048 


0-827 


1-0 


-3365 


0-507 



multiplying Column II of the Table on p. 301 by 1*3125, which is now the 
proper ratio required (instead of 1*348) to convert the tilting effects there 
considered into total deflections. This result might have been foreseen. 
Eecurring to the typical formulae (91) and (92) we note that 



2tt 7 H 



gU y 4:iTyfM 



#H aft 



9 2/* sr 
Hence we find after a few reductions 



> («-/3) 5 



(105) 



dx\ 



&' 



/3 



gj (*-/3f 2/* L a 3 + /S 3 m + a * 3 + /3 3 m-/3. 



sm mx. 



(106) 

If /3 = \a the last term in the bracket is small, and the expression is almost 
in a constant ratio to that which occurs in (50) provided tc be written for a. 

The deflections at given distances are, however, diminished from the values 
based on Art. 5 above, where undisturbed gravity was alone regarded in 
calculating the strains f and at great distances the effect is considerable. 
Since on the present reckoning, 

. a-/3_4 



K 



a 



a 



5 



any line in the above Table relates to a distance four-fifths as great as that 
implied in the corresponding line of the Table on p. 301. To take a 
numerical example, let k, x = 0*8, which means a distance of about 2300 km. 
The total deflection is then only 0*507 of that given by a purely elastic 
theory ; on the basis of Art. 5 the ratio would be 0*579. 

8. The investigation for the case of a load M concentrated at a point of 
the surface follows a parallel course. We should find 

Mg f°°m 2 — a/3m/(q— /9) 



w = 



J (mr) dm 



4 iTfju J o (m -f a) (m — ft) 



9 



M7 ri 

g \r 



a + PJo m+* <* + #. 



00 



Jo (tnr) 
m—ft 



dm 



(108) 



Vertical by Tidal Loading of the Earth's Surface. 309 

The principal value of the indeterminate integral is, as shown in the 
Appendix, 

— ^ — zfdm = 7rK (/3r)— ■ — ^—~ dm 

^ J Jo m + fi 



m—ft 



\ ttKo C8r)- i 0r-&& + p^p - • • ■ h • (109) 



From (107) we derive 



a" 



<fa>o = Mff /l 

<#r 4 7r^ \r 2 a 2 — /3 2 



f — ^dm — —^ V^n ( 110 ) 

o ffl + « a -p Jo m ~ P J 



with a corresponding formula for V . The former of these integrals is 
evaluated in (74). As regards the latter we have by differentiation of (109) 



1 t-. f °° mJi (mr) 7 
Jo m—/3 



i7rK 3 (/3r)+i 1-^? + (/3r)4 



l 2 . 3 ' l 2 . 3 2 . 5 



(111) 



j 



The following Table, calculated on the numerical assumptions adopted in 
the preceding section, gives in Column I the ratio which the surface 
depression w bears to the value given by the first term in (107). Column II 
gives the ratio of the tilting effect to that given by the first term of (110), 
whilst Column III gives the ratio of the total deflection to the same term. 
It is assumed that the use of the " principal values " of the indeterminate 
integrals can be justified by considerations similar to those adduced in the 
two-dimensional case. 



ar. 


I. 


II. 


III. 


ar. 


I. 


II. 


III. 


o-o 


1-0 


1-0 


1 -3125 i 


0*5 


-3735 


-6334 


0*860 


0-05 


*8339 


'9499 


1 *246 | 


0'6 


-3265 


*5840 


0-799 


o-i 


-1352 


-9039 


1 -193 ! 


0-7 


'2871 


-5392 


0-744 


0'2 


-5992 


-8219 


1 '092 


0-8 


-2542 


-4996 


0-696 


0*3 


'5037 


-7508 


1 -005 ; 


0-9 


-2260 


-4627 


0*651 


0*4 


0*4312 


-6886 


-928 


1*0 


-2020 


-4297 


0-610 



The comments made at the end of the preceding section apply equally 
well to these results. Thus, if kt = 0*8, the total deflection is 0*610 of the 
tilting effect given by a purely elastic theory, whereas on the assumptions of 
Art. 6 the ratio would be 0*690. 



Appendix. 

The following are proofs of some of the formulae employed in the paper. 
I do not expect that they are altogether new, but the literature of BesseFs 
functions (for instance) is so diffuse, and the proofs so entangled, that it is 
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sometimes difficult to find a concise demonstration, even of an important 
formula. 

(A) The equation 

f'(v) + lf(xj+f(v) = l (112) 

is satisfied by 

and also by 



•00 



f(x) = \ e- xahhw du, (114) 

o 



x being of course assumed to be positive. For (113) makes 

/'(^)+^/'(^)v(^) = } o w { ( ^-^^+^}jo(.)& 

' o W {j "(,)-ij„'(.) + Jo(.)}^ 

1 TOO I 

= _± J '(z)dz = -, (115) 

by partial integrations, and in virtue of Bessel's equation. Again (114) 
makes 



/'» +;/'(*) +/(«) = 

yV 



g -.r smh m COiSn 2 ^ _ _ g^fr w ^ 
\ X J 

If 00 1 

= -- ^ (<r*atahu cosh ?/,) = '-. (116) 

a? J o ^ 

The difference between the integrals in (113) and (114) must therefore be of 

the form 

AJo(#) + BKo(;c), 

where K is defined by (73) ante. Now the leading terms in the asymptotic 
expansions of J (x) and K (x) for large values of x are*. 



\2 



( — V sin (a; + 1 7r) and ( — )" cos (a? + \ 7r), 

\7raV v 4 7 \7raj/ 

respectively, whereas either definite integral is evidently at most of the order 
1/x. It follows that A = 0, B = 0, and therefore 

r hM dz = rg-xsinhu^ (H7) 

Jo 2 + ^ Jo 

as stated in (71). 

* c Hydrodynamics,' Art. 194. 
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(B) Consider the integral 

where f= g + vrj. Let this be taken round the contour of the infinite region 
lying above the axis of £, the point f= — ft being excluded by a small 
semicircle described about it. The infinitely distant parts of the contour 
contribute nothing to the final result, and we find 

P 



_»£+£ 



df _,rie-*«+ f° pf^r d£ = 0, (119) 

JO £~TH 



(*oo p—i^x fee pi'"x 

^o^Ho^- 7 ^ (120) 
Hence, separating the real and imaginary parts, 

COS cJu t c » r\ COS QtAy -i t,' , _, ~. _, . 

fe d£ = — -n- sin /3x + fc—o^' (121) 



P 



o Z-0 b ^ Jo f + /8 

J f^f ^ = t cos fix- |" ^| rff . (122) 

The latter formula has been assumed in (109). 

(C) If in (121) we write rcos0 for x, and integrate with respect to 6 
between the limits and J 7r, we obtain 

p f 00 MjQ^ = _ 2 [ h \ in (i3 rQO a0)de+ r^fW-d^ "(123) 
Jo g — p Jo ' Jo g-rp 

2 r^ 71 * 

since Jo (f;r) = - cos (f r cos (9) <tf#. (124) 

Again, writing r cosh^ for a?, in (122), and integrating from u = to 
^ = oo , we have 

p | "f=§^ = 7rK »^')-|r^9^ a 25 ) 

in virtue of (73), and of Mehler's formula 

2 f 00 
Jo(£r) = — ! sin (£r cosh u) du. (126) 

^Jo 

By combining the formulae (123) and (125), we obtain another proof 
of (74). 

(D) To sum the series in (102) we take the integral 

cos (| & — x) £ 



i 



/*, - i r d t ( 127 ) 



round the contour of the infinite region lying to the right of the axis of r), 
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the poles of the integrand being excluded as usual by small circles drawn 
round them. These poles are at the points 

f = sir J a, [s = 1, 3, 5, ...], (128) 

and the value of the integral taken in the negative direction round any one 

of them is 

4:7ri sin (s7rx/a) ^ ~q\ 

a S7r/a-\-cc 

If x be positive, and <%a 9 as we will suppose, the part of the integral due 
to the infinitely distant portions of the boundary will vanish. The sum of the 
series of terms of the type (129) is therefore equal to 

f°° cosh (ha—x)7) : . 7 /-toA\ 

,. , \ 2 — 7rT ± -%cbj. (130) 

J_oo W + a) COSh f at] K 

This leads to the formula (102). 
In the case of the integral 






there is an additional pole £ = /3, which gives a term 



2nri 2£i£M=£), (132) 

eos|/3& 



The formula (103) easily follows. 



